A single structured light beam as an atomic cloud splitter 
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Abstract 

We propose a scheme to split a cloud of cold non-interacting neutral atoms based on their dipole 

O 

interaction with a single structured light beam which exhibits parabolic cylindrical symmetry. 

^ Using semiclassical numerical simulations, we establish a direct relationship between the general 

1—5 

ly-^ properties of the light beam and the relevant geometric and kinematic properties acquired by the 

atomic cloud as its passes through the beam. 
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An optimal implementation of an atomic thermal cloud splitter requires achievement of 
a predetermined efficiency and deviation angle, as well as prevention of heating processes. 
Optical devices for splitting a free falling atomic cloud usually rely on two, or more, non- 
collinear [Ij or collinear [2j light beams, since a single Gaussian laser beam orthogonal to the 
falling atoms yields deviations of only some ten milliradians [3]. On chip, optical splitting 
requires either a pulse sequence of tightly hnearly polarized beams [4J or the splitting into 
two of a former one well dipole trapping beam via acousto-optical modulation [5]. It is also 
possible to use on chip magnetic fields to split the atomic cloud by including current carrier 
wires with varying bias or distance between them [21 [7] . 

Here, we propose to reconsider the simple configuration of a single but structured laser 
beam orthogonal to a free falling atomic thermal cloud. Light beams with special intensity 
or phase structure yielding peculiar dynamical properties have been extensively used for the 
manipulation of microparticles [8] and cold atoms [9]. Most theoretical and experimental 
studies for cold atoms focus on using such light beams as wave guides for atom transport 
with the possibility of exchanging orbital angular momentum P [TOl [HI [121 [131 [IH [13, 
[16] or other dynamical properties [1^ between the light and the atoms. We shall show 
that experimentally accessible structured light beams can also be used efficiently as atomic 
thermal cloud splitters. This proposal is based on semiclassical calculations that allow an 
estimate of the relevant geometrical and mechanical properties of the atomic cloud, such 
as deviation angles and kinetic energy, in terms of the atomic and light beams setup. In 
these calculations, the center of mass motion of each atom is described using Newtonian 
mechanics with the optical force derived from quantum theory |27j . 

The structured light beams we study are known as vector Weber beams. Under ideal 
conditions, the intensity pattern of an electromagnetic (EM) Weber beam is invariant along 
the main direction of propagation. Its transverse structure is naturally described in terms 
of parabolic coordinates u and v related to the Cartesian coordinates through the equations 
X = {u^ — f^)/2, y = uv, ue{—oo,oo), ve[0,oo). Its associated electric and magnetic fields 
can be written as a superposition of transverse electric (TE) and magnetic (TM) modes [18] 

= d,t (^(^^)N - ^(^^^)M) m^. (1) 
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The vector operators M and N are 

M = h'^dct {eudy - e^du) , 

N = h-^d, {eudu + e,d,) - e,Vi, (2) 

where h = \/v? + v'^ denotes the scahng factor, e^, C, = u, v, z, the unitary vectors cor- 
responding to parabohc-cyhndrical coordinates, and the shorthand notation dx is used for 
partial derivatives with respect to the variable x. The scalar field is the solution of the 
wave equation: 

= f^p,fcx,a(M)Vp,fc^,a(^^)e V J ^ (3) 



where 



2 "p-l frin a Tin J o 
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with Up = 1 {up = 3 ) for even (odd) parity functions, iFi{xi,X2,X3) is the confluent hyper- 
geometric function directly related to Weber functions fl9[ |20| [21]. Up to a normalization 
factor, the scalar wave function 'ipp^k_i_,a in terms of the Weber beam angular spectra 21 is 
given by [22]: 



gia In I tan ip/2\ 

2le(a;<^) = (5) 
^, , . J i^e{a;^) <^e(-7r,0) 

2lo(a;(/?) = < (6) 

The labels set k = {p,kz,uj,a} in Eq. ([T]) refers to the mathematical parity, p even or odd, 
2;-component of the wave vector, kz, frequency, uj = ck = c^Jkl + k\, and continuous order, 
a G M. The Weber beam transverse structure exhibits extended zones where the field has 
zero, or close to zero, values as illustrated in Fig. ([T]). For the order parameter a 7^ 0, the low 
intensity zone is to the right (left) depending on the parameter a being negative (positive); 
its area increases in size as the absolute value of a increases in magnitude. This dark region 
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is approximately delimited by the parabola y = a/2m|^(|x| — where um is the value 

of u at the first maxima of the function ?7p,fcj^,a(w). In Fig. ([2]), um as a function of a is 
illustrated. 

The dynamical constants of Weber photons are [23]: its momentum along z, hkz, its 
energy hu, and its eigenvalue, h'^k±a, for the quantum operator 

^ " {IzPy + Pyl 

1=1 



(7) 



Here N^f^ denotes the number operator of the k mode, Iz = —i{fx V)^ the generator of 
rotations along the z-axis and Py = —idy the generator of translations along the |/-axis. Then, 
it would be said that the parameter a determines the product of the angular momentum 
along z and the momentum along ?/ of a Weber photon whenever it could be shown that 
charged particles exchange that mechanical property with those photons. 
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Figure 1: (Color Online) Sampler of the transverse structure of the intensity of Weber EM fields 
as defined by Eqs. ([T]). Even/odd EM Weber fields with eigenvalues (a/d) a = —2, (b/e) a = 0, 
(c/f) a = 5. 
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Figure 2: (Color Online) Value for um, the coordinate that defines the parabola delimiting the low 
intensity region of a Weber field, as a function of the eigenvalue a. The length unit was taken as 
the wavelength of the beam and kz = 0.995fc. 

Recently, Weber beams of zero-order were experimentally generated by means of a thin 
annular slit modulated by the proper angular spectra [21]. This setup was conceived as a 
variation of that originally used by Durnin et al.[25] for generating Bessel beams. Higher 
order Weber beams can be produced using holograms encoded either on plates |21] or on 
spatial light modulators [SB] . 

The interaction between a two level atom and an EM Weber wave with a frequency close 
to resonance is determined by the coupling factor 

where the ± super index refers to the two possible cases for the atomic dipole transition 
element /xf2 = l^t^vii^x^f'^y)- In a semiclassical treatment, the gradient of 5^, V(7 = {oL+il3)g^ 
defines the force experienced by the atom. The expression for the average semiclassical 
velocity dependent force valid for both propagating and standing beams is |27j : 



hV 



[{D{1 - p)v ■ (x) + r/2]/3 + [{v ■ (3) - 5uj]cx 



with r the Einstein coefficient, F = 4A;^|/xi2p/3/i, 5uj the detuning between the wave fre- 
quency uj and the transition frequency uiq, 5uj = uj — Uq, p = 2|(7p/ ((r/2)^ + (5a;^) a saturation 



parameter linked to the difference D between the populations of the atom two levels, D = 
l/(l+p), and finally = 2|^|V|7f, with 7' = {v-a){l-p){l+p)-^ + T/2 + t{-5uj + {v-l3)). 

Our proposed scheme is focused on the red detuned far off resonance case so that non- 
conservative terms arising from the velocity dependence of the force are not dominant [28] . 
Nevertheless, we keep the velocity dependent terms in the numerical calculations in order 
to prevent disregarding potentially relevant effects, since the atom may increase its kinetic 
energy as it interacts with the light beam. The numerical simulations consider a TE laser 
beam detuned 67 nm to the red of the 5 ^6*1/2 - 5 transition at 795 nm of ^^Rb with 

irradiance in the range ~ 1.5 — 6 W/ cm^. As natural length and time units, we take the 
laser wavelength and the inverse of the Einstein coefficient F, which is 3.7 x 10^s~^ for the 
state 5 2Pi/2 of ^^Rb. 

An atomic cloud is released into a region where a Weber beam propagating in the horizon- 
tal direction was generated. The atoms have an initial velocity oriented by the gravitational 
field in the negative x-direction, with a random smaller component in the y or z-direction 
(we shall take \vx\ > 10v± = lO{Vy + v'^Y^'^). The Weber parameter a is taken negative so 
that, initially, the atomic cloud faces the low intensity region; for x > 0, due to the red 
detuning of the light, the atoms that arrive in the dark region of the beam are attracted 
towards the delimiting parabola mentioned above and the splitting of the cloud takes place. 
In the region a; < these atoms find a series of effective potential well channels that give 
rise to their focusing towards those channels. The atoms that arrive to the Weber beam 
outside the central dark region experience either a strong deflection when their initial po- 
sition is on the bright parabolas, or a weak deflection if their initial position is outside 
those parabolas. This behavior was found to be generic in all numerical simulations we 
performed. For the parameters mentioned above and mv1/2kB ~ l/^K, the irradiance can 
be chosen so that the maximum deflection angle for the atoms that arrive in the dark region 
is 0^^^ ~ arctan((u|_^/|a;o|(l — u\j/2\xq\)Y^'^)^ where Xq is the initial atom x-coordinate and 
um is the coordinate that defines the parabola delimiting the dark region; um is determined 
by the a parameter of the Weber beam. Fig. ([2]). As expected, higher (lower) initial atomic 
kinetic energies require higher (lower) irradiance to achieve such deflection angles. 

These results are illustrated in Fig. ([s]) for a = —5 and Xq = 150A. If mvl/2kB ~ 1.5/iK 
and the irradiance is 1.725 W/cm'^, the atoms that arrive to the beam within a distance of ~ A 
to the X-axis are violently deflected, while those arriving in a zone X < \y\ < 2u\,^ ~ 50 A have 
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Figure 3: (Color Online) Evolution of an atomic cloud as it interacts with an odd TE Weber beam 
whose intensity pattern is illustrated in black and white. The atoms initial velocities are mainly 
in the direction of the gravitational field, Vx ~ —0.6 mXT with random Vy and Vz components, 
l^^xl ^ 10f_L, corresponding to T ~ 1.5 fiK. The Weber beam is characterized by a = —5, 
kz = 0.995 k and irradiance (a,b) 1.725 W/cm? and (c,d) 6 W/cm^. The atoms motion takes place 
within a ~ 20 A wide region in the z-axis. 

deflection angles up to 6'^"^ ~0.17 radians, Fig. (|3^); those arriving at 80A < \y\ < 150A 
move almost vertically. Fig. ([sjo). If the irradiance is almost quadrupled the deflection 
angle reach values of the order 0^°-^ ~0.35 radians for atoms arriving at the dark zone, 
Fig. while the atoms arriving at 80A < \y\ < 120A are partially focused towards the 
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X-axis, Fig. (|3]i). This effect is also interesting by itself since it could be used for merging 
independent clouds. In Fig. (|4^), the deflection angle, 9d, is illustrated for all arriving zones. 
Typical values of 6^ are tenths of radians for atoms arriving at the dark zone of the beam. 

Due to the parabolic symmetry of the problem, the product of the atomic angular momen- 
tum Lz and the linear momentum component Py, Aatomic = LzPy, is a relevant dynamical 
property for the description of the atomic cloud. As mentioned above, Weber beams carry 
a well defined value of the similar EM property Aem, Eq. ([T]). Through its classical motion, 
each atom locally detects the density dAEAi/dV, and could modify its Aatomic value. For 
the examples illustrated in Fig. (|3^,b), initially the atoms have Aatomid^) ^ 150h?k± and 
finally highly deflected atoms have 10^h^k± < Aatomic < lO^fi^kx as shown in Fig. ^p). A 
detailed analysis lets us see that atoms arriving outside the dark beam region may exhibit 
high fluctuations in their Aatomic value as a function of time, while atoms focused into the 
same bright parabola for x < have similar asymptotic Aatomic values. These results show 
that Weber photons may transfer the dynamical variable A to atoms. 

As for the kinetic energy acquired by the atoms, we divide it in two parts: the contribution 
due to the velocity in the gravitational direction Kx = and that in the transverse 

direction Kyz = mv\/2. For the example illustrated in Fig. ([5^), initially Kx{0)/kB ~ l.SyuK 
and Kyz{0)/kB ~ 0.02yuK finally K^/kB ~ 20/iK and Kyz/kB{0) < 1/iK for atoms arriving 
at the dark region, and Kyz{0)/kB ~ IfiK for atoms arriving in the other regions. The 
movement in the x-direction is mainly driven by the gravitational field. The spreading 
in the kinetic energies is smaller for atoms arriving at the dark zone. At short times, 
< t < 1.5 X lO^r, \vx\ diminishes while in the average v± ~ \vy\ increases. That is, the 
atomic deflection takes place mainly in that time interval. Thus, if smaller final kinetic 
energy for the atoms is desired, the region of the optical field could be shortened along the 
X-axis once the atomic thermal cloud has been split. 

Summarizing, we have presented a scheme to split a cold non-interacting atomic cloud 
using a single Weber beam. This scheme requires no further optics or electronics but that 
used to generate the Weber beam. The atomic deflection angles are a function of the 
Weber parameter a and the irradiance of the beam, as well as the initial kinetic energy and 
arriving zone of the atoms. We have shown numerically that it is possible to manipulate 
the population of the split components of the atomic cloud through the control of the entry 
point of the cloud with respect to the center of symmetry of the EM field intensity pattern. 
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Figure 4: (Color Online) (a) Absolute value of the deviation angle, Q = arctany/x and (b) temporal 
evolution for the values of the mechanical property Anatomic = LzPy for different components of the 
clouds presented in Figs, (pk-b) 



In fact, we have shown that this arrival dependence may be used to convert the hght beam 
into a merging device. The atoms, when guided by a parabohc bright region of the beam, 
acquire values of Aatomic = L^Py with small dispersion. The main source of kinetic energy 
for the atoms is provided by the gravitational field. 

For atoms that initially have kinetic energies in the quantum regime, lower irradiance of 
the light beam is expected to achieve similar deflection angles as those reported here. Under 
such conditions, studies on the quantum evolution in phase space of the atomic beam and 
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Figure 5: (Color Online) Contribution to the kinetic energy from, (a) velocities parallel to the 
gravitational force direction, Kx/ks = mv^/^kB, and (b) velocities perpendicular to the gravita- 
tional force direction, Kyz/ks = mv'j_/2kB, for different components of the clouds presented in 
Figs, ^-h) 

its consequences on interferometric experiments would be especially interesting due to the 
topological structure of the EM Weber field |23] ■ 
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